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Abstract The energy levels of quantum systems are determined by quantization conditions. 
For one-dimensional anharmonic oscillators, one can transform the Schrodinger equation into 
a Riccati form, i.e., in terms of the logarithmic derivative of the wave function. A perturbative 
expansion of the logarithmic derivative of the wave function can easily be obtained. The Bohr- 
Sommerfeld quantization condition can be expressed in terms of a contour integral around the 
poles of the logarithmic derivative. Its functional form is B m (E,g) = n + |, where B is a 
characteristic function of the anharmonic oscillator of degree m, E is the resonance energy, and 
g is the coupling constant. A recursive scheme can be devised which facilitates the evaluation 
of higher-order Wentzel-Kramers-Brioullin (WKB) approximants. The WKB expansion of the 
logarithmic derivative of the wave function has a cut in the tunneling region. The contour 
integral about the tunneling region yields the instanton action plus corrections, summarized 
in a second characteristic function A m (E,g). The evaluation of A m (E,g) by the method of 
asymptotic matching is discussed for the case of the cubic oscillator of degree m = 3. 
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1 Introduction 

The cubic anharmonic oscillator, defined by the Hamiltonian, 



is a paradigmatic example of a quantum mechanical problem which gives rise to complex resonance energies. A 
particle initially trapped in the region q m may tunnel through the classically forbidden well and escape toward 
q — > —oo. For positive coupling g > 0, the energies can be determined numerically by the method of complex 
scaling[2J [23]. One scales the coordinate as q — s- qe l6 for the cubic oscillator, which results in the Hamiltonian 



h ^ = -2W + 2 q2 + ^ 9q3 



(1) 




(2) 
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If one diagonalizes this complex operator in the basis of harmonic oscillator wavefunctions {(p n (q)}n=o x • f° r l ar 9 e 
enough N max , one can numerically determine the (complex) resonance energies of the original cubic Hamilto- 
nian fll]). As discussed in Refs. [TJ [24], these resonance energies are independent of 9, provided we choose 9 
sufficiently large so that the rotated branch passes the position of the resonance under investigation. The numeri- 
cal variation of the resonance energies of the cubic oscillator, 

e£%) = Rc (e^(g)) + ilm (e^ (g)) , (3) 

as a function of g under a suitable increase of N max can be used to estimate the numerical uncertainty of the 
numerical results for the complex resonance energies [16j . Throughout this article, we keep the superscript "3" of 
the resonance energy eh\g) of the cubic oscillator in order to remain consistent with the notation of a brief Letter 
(see Ref. [18]), in which part of the results presented here have already been indicated. 

For reasons which will explained below, we associate the resonance energy (with a negative imaginary part) with 
a value of the coupling parameter g that has an infinitesimal positive imaginary part, whereas a coupling with an 
infinitesimal negative imaginary part is associated with an antiresonance, that has a positive imaginary part. The 
resonance energies of the cubic potential as a function of g have a branch cut along the positive real axis, and 
fulfill the dispersion relation [HE], 

2 n J s (s - g) 

We note that the spectrum of resonances is invariant under the transformation y/g — > —\/g in Eq. ([1]), which can 
be compensated by a parity transformation q — > —q (the latter leaves the spectrum manifestly invariant). 

For weak positive coupling g = \g\ — > (with Img > 0), the imaginary part of the resonance energies is determined 
by a nonanalytic factor, which reads [16] 



Im4 3) (.g + i0) 2 "g-n-i/* e ^(JL\ for g = \g\^0 t lm.g>0. (5) 

n'.y/Tr \ log/ 

For large positive coupling g = \g\ — >• oo + iO, the first three resonance energies possess the following asymptotics 
[see Eq. (17) and Table II of Ref. [16]], 

4 3) (<? + i0) - g 1/5 (0.617 160 050- i 0.448 393 023) for g = \g\ oo + iO , 
4 3) (c/ + i0) - g 1/5 (2.193 309 731 -il. 593 532 797) for g = \g\ -> oo + iO , 
4 3) ( 5 + i0) ~.g 1/5 (4.036 380 020 -i 2.932 601 744) for g = \g\ -> oo + iO . (6) 

All of the above results have a complex argument of — n/5, as they should. For intermediate g, a large number of 
resonance energies can be determined numerically by a suitable basis-set method and be used in order to construct 
a complex-scaled time propagation algorithm for a quantum mechanical wave packet (see Sec. Ill of Ref. |16|). 

A second available method for the calculation of resonance energies consists in the resummation of the divergent 
perturbation theory in complex directions of the parameters. The perturbative coefficients e^' K in the perturbation 
series 

oo 

$H9)~J2£!k9 K , 9^0, (7) 

K=0 

which describe the nth resonance energy of the cubic oscillator, diverge factorially for large perturbation theory 
order K. Indeed, for large K, the leading asymptotics read (see, e.g., Ref. 1181), 



r 



(3) 

n,K 



22n-l/2-K ^n+K+1/2 

it 3 / 2 nl 



r (n + K + ±) , K -> oo . (8) 



With the exception of the coefficient of order zero, which reads e^ 3 Q = n+1/2, all e^ K are negative for K > 1. For 
g > 0, the series ([7]) is formally not Borel summable. However, numerical evidence (TO] [T2j [15] [16] suggests that 
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the Borel-Pade summation method with the Laplace-Borel integration being carried out along complex contours, 
as given by Eqs. (196)-(198) of Ref. \W\, provides useful numerical approximants to the resonance energies. 
The method has been put on rigorous mathematical grounds recently \6\ in the framework of distributional Borel 
summability [8l[9]. When the integration contour C+i in the conventions of Ref. [10] is employed, a negative 
sign is obtained for the imaginary part of the resonance energy. The accuracy obtainable using complex Borel 
summation of the weak-coupling expansion is restricted by numerical oscillations in the values of the complex 
Borel-Pade transforms as the order of the Pade transformation is increased. These oscillations cannot be overcome 
even if extended-precision arithmetic is used in intermediate steps of the calculation. These oscillations have been 
described in Refs. [161 [21] f° r the complex Borel-Pade resummation method. Indeed, it has been observed in [16] 
that, at a relatively moderate coupling g = 0.6, the ground-state energy of the cubic oscillator as determined by 
resummation cannot be calculated to better accuracy than 



by Borel resummation. The oscillations, which are described by the numerical uncertainty in the last digit reported 
in result of Eq. ©, cannot be overcome when the (Borel) transformation order is increased and appear to represent 
a fundamental limit of the convergence of resummed weak-coupling perturbation theory in the case of a moderate 
and large (modulus of the) coupling parameter g. A more accurate result obtained by complex scaling for the 
same coupling parameter is 



The applicability of the weak-coupling expansion for the cubic oscillator resonance energies thus is found to be 
severely limited in the domain of moderate and large coupling, even if the perturbative expansion is a posteriori 
enhanced by a summation method acting in the complex plane. Likewise, the leading-order analytic result ([5]) 
cannot be applied to the domain of moderate and large coupling g. In this article, we explore the question whether 
one can find higher-order analytic formulas for the corrections to the leading result ([5|) for the decay width of the 
ground-state resonance energy (and, possibly, an arbitrary excited resonance energy) of the cubic oscillator, and if 
yes, how large the correction terms are. 

Indeed, in Ref. [18], an analysis has been presented which allows us to calculate higher-order correction terms 
to the resonance energies for even and odd anharmonic oscillators of arbitrary degree, and to formulate modified 
Bohr-Sommerfeld quantization conditions for anharmonic oscillators such as the cubic one. Here, we provide a few 
more details of the analysis that leads to the results of Ref. [18] , with a special emphasis on the cubic potential. To 
this end, we review in Sec.[2]the calculation of higher-order perturbative (in g) approximants to the wave function 
for the cubic potential. In Sec. [3] the calculation of higher-order corrections to the WKB approximation of the 
wave function is reviewed. The modified quantization condition pertaining to the cubic oscillator is described in 
Sec. [4] with a few new results for higher-order terms reserved for Sec. [5] Finally, conclusions are drawn in Sec. [6] 

2 Perturbative expansion 

The Schrodinger equation corresponding to the Hamiltonian given in Eq. ([1]) reads 



e£>(g = 0.6) = 0.554(1) - 0.351(6) i 



(9) 



e£> (g = 0.6) = 0.554 053 519 - 0.351 401 778 i . 



(10) 




(11) 



If we transform the coordinate according to q — >• q/^/g, we obtain 





(12) 



We may transform to the Riccati equation by setting 



V"{q) S*(q) s'(q) 



(13) 



9 



f(q) g 2 g 
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The Riccati form of the Schrodinger equation then reads, 



gs'(q)-s 2 (q)+U 2 (q) = 0, U(q) = y/2[V(q) - gE] . (14) 

The zeroth-order term is 

«(g)«*o(<z)=«(«) = V2Vfa). (15) 
The following recursion yields higher-order approximants (in g) for the logarithmic derivative of the wave function, 

s(?) = Yl 9 K s K (q) , s (q) = u(q), si(q) = - > 

We can divide s(g) into a symmetric and an antisymmetric component under the simultaneous interchange g — > —g, 

E -> 

s (<7) = M?) + M<7) = s+{q,E,g) + S-(q,E,g) , s±(q,-E,-g) = ±s ± (q,E,g) . (17) 
From the following system of equations, 

gs'_(q) - s+(q) 2 - s_(g) 2 +U 2 (q) = , (18a) 
5 sV(g)-2 S+ (g) S _( g )=0, (18b) 
we can eliminate the antisymmetric component by setting 

•-(*) = | (19) 
The Bohr-Sommerfeld quantization condition reads 

-2k-gl dZS+{z)=n+ l (20) 

where C is a contour that encloses the zeros of the function s+(z) in the anticlockwise (mathematically positive) 
direction, and n is the quantum number of the level. The function B 3 (E,g), defined by 

B 3 (E,g) = -^— £dzs + (z), (21) 

[we recall the implicit dependence of s + (z) on g and E\ provides for a universal means of determining the perturbative 
expansion for an arbitrary excited level of the cubic potential, by means of the perturbative quantization condition 
B 3 (E, G) = n + |. A calculation leads to the result 

7 , 15 t^A , i( 1365 ™ , 1155 r&\ , 3 / H9119 , 285285 ^ , 255255 



B 3 (E, g) = E + g — + — E 2 +9 S + E 3 ) + a 3 [ + + + C(g 4 ). 

a| 16 4 J y V 64 16 / y V 2048 256 128 / w ; 

(22) 



3 WKB expansion and contour integral 

The WKB expansion is an expansion in g at g E fixed, i.e., it implies large values of E. We use this expansion in 
order to derive, roughly speaking, higher-order corrections to the instanton action which can later be used in order 
to derive higher-order corrections to resonance energies E of levels whose principal quantum number is not large 
(this is based on a rather fortunate enhancement of the applicability regions of certain expansions, as explained 



4 



below). The instanton configuration for the cubic potential has been discussed in Ref. [17], with a graphical 
representation being given in Fig. 3(b) of [17]. In order to achieve an instanton configuration pertaining to positive 
values of the coordinates, we scale the coordinate as q — > —q/^fg. We start from 

\l® + \* + °t)+ = E + (23) 
and transform according to q — > —q/^fg, which gives 

+ W(q)) <t> = gEcf>, W{q)= l -q 2 - q 3 . (24) 



(25) 



2 dq 2 ' ""'J r a ~~ r ' " w 2 
We note that W(q) differs from V(q) in the sign of the cubic term. Setting 

4/ _ S(g) 

where the logarithmic derivative is denoted S(g) instead of s(q), we obtain 



5 S'(g)-S 2 ( g )+T 2 (g)=0, T{q) = ^2[W{q)-gE]. (26) 
The zeroth-order term is different as compared to the perturbative expansion, 



S(«) «S (g) = V2[W(g) -<?£]. (27) 

A recursion in g can be defined, 



S(g)= X>*S*(g), So(g)=r(g), Si(g) 



k:=o 



2S (q) ' 



Like the perturbative expansion, the WKB expansion of the logarithmic derivative of the wave function can be 
calculated recursively in ascending powers of g, and we can also distinguish between even and odd components 
under the symmetry operator g — > — g, E — > —E, 

S(q) = S+(q) + S_(g) = S+(q, E, g) + S_(g, E,g) , S±(g, -<?) = ± S±(g, £, «?) . (29) 

The contour integral about the cut of the function S + (g) leads to a more complex structure as compared to the 
right-hand side of Eq. (f20j) . It reads 

(30) 



i fdzS + (q)=A 3 (E,g) + ±ln(2Tr)-]n[r(±-B 3 (E,g))] +B 3 (E,g) ln(-|) . 



Here, C is a contour that encloses the cut of the WKB expansion of the logarithmic derivative of the wave function 
in the clockwise sense. We thus have to calculate the integral 



T = - g j> dq S+ (q) « i y dq [S (q) + g 2 S 2 (g) + <? 4 S 4 (g)] = T + T 2 + T 4 



(31) 



Here, the WKB approximants So(q), S 2 (g) and 84(g) are defined in Eq. (|28|) , and the integrals T 0l T 2 and T4 
correspond to the contour integrals of the WKB approximants of the respective order. The zeroth-order term is 

T = -IdqS {q) (32) 
c 



5 



which involves the zeroth-order WKB approximant 



S (g) - y/2W(q) - gE = vV - 2g3 - gE . (33) 

If we define the square root function to have its branch cut along the positive real axis, then the above expression 
has its branch cut (in the limit gE — > 0) from q = to q = |. Directly above the cut, the value of the square root 
is positive, while directly below, it is negative. As the contour encircles the cut in the clockwise direction, we have 
an integration interval from zero to i above the real axis, while below the real axis, we go from | to zero. With 
this definition, the contour integral is 

1/2 

To = - j> dq vV - 2q 3 - gE = - Re J dq (q 2 - 2q 3 - gE) 1 ' 2 . (34) 
C 

For g E > 0, the cut of the WKB approximant does not extend along the full interval (0, \) any more. However, 
by specifying the real part of the final result, we pick up only those terms which are due to the cut. 

The method of asymptotic matching (see Ch. 7.4 of Ref. [4]) can be used in order to evaluate (j34|) . It involves an 
overlapping parameter e, that is also used in Lamb shift calculations [TTJ[T4j[20] in order to separate the low-energy 
from the high-energy contribution to the bound-electron self-energy. The overlapping parameter fulfills 



0<9£«e. 



(35) 



We then separate the integration interval into three parts, the first of which extends from (0,e), the second, from 
(e, | — e), and the third complementing interval finally reads (i — e i). 



The method of evaluation depends on the particular integration interval. For the interval (0,e), we have to be 
careful to avoid divergences for very small q = 0. We have to keep the expression y/q 2 — 2g E within So(q) in 
unexpanded form. However, we can expand the expression \J 'q 2 (1 



2q) - 2gE = y/q 2 -2gE-2q 3 in the q 3 
term, then do the q integration, subsequently expand the result in gE, and finally in e up to order e°, keeping all 
divergent terms [11] [14j [20] . This leads to the result 



h=Ref dq (q 2 - 2q 3 - gE) 1 / 2 = (gE) 



2f_ 

gE 



(gE) 2 



47 
16 



1 

4^2" 



27 



15 



In 



— ^ 

gE) 



(gE) 3 
(gE) 4 



13327 1 5 35 105 1155 j/26 2 
192 + 8^ + 6^ + 8^ + ~e 32~ " \g~E 

6364777 5 7 315 385 15015 
3072 + 48^ + 8^5 + 16? + 16? + 128e 2 4 



45045 255255 



64e 



256 



gE 



0(gEf 



(36) 



For the interval (e, i — e), we can expand in gE to any power, because all divergences near q = or q = 1/2 are 
cut off by the overlapping e parameter. The result reads 



1/2-6 



h = Rc 



dq(q 2 



gE) 1 ' 2 



+ (gE) 3 

_7_ 

~8? 



23189 
192 

315 



1 

8? 



_5_ 

6e» 



1 

15 " 

35 
8? 



(9E) 
105 

" 17 



In 



4^2 



- (gE) 2 

40^2 1155 



47 



385 15015 45045 



3e 3/2 £ l/2 

2\/2 U0V2 



64e 4 16e 3 128e 2 



64e 



F 5/2 



F 3/2 



32 
1120V^ 

f l/2 



1 

4? 



3 



2y/2 15 



In 



In — 



4 

255255 
256 



(gE) 4 



( - 



5241655 



1536 



48e 6 

0(gEf. 

(37) 
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The result from the interval I 3 compensates the divergences of fractional order in e from I 2 , 

1/2 

1. /1 

(9E) 3 

c-/ - 

1/2-e 

' 2\/2 140V2 1120^2 



/ 3 = Rc y d<? (g 2 - 2g 3 - s^) 1 / 2 = 

C. 



4^2 40^/2 

3^2 + 7^" 



£ 5/2 3e 3/2 



a/2 



G(gEY 



In the result for T , 



T = - (h + I 2 + I 3 ) 
9 



the overlapping parameter e cancels, and we obtain 
2 



T ° 15g 
+ 9 2 

The second term reads 



3043 1155 - , 
E 3 H £ 3 In 



141 9 15 n 
£ 2 E In 

8 4 



1G 







(f): 


+ .9 3 



5616029, 255255 . ( gE 
E H E In 



512 



128 



V 8 



T 2 = dqS 2 (q) = 



16431 „ 1365 
E 



128 



64 



41 


7 




16 


(¥ 


)] 







In 



862501 
128 



-£ 2 



285285 
256 



E< In [ ^ 



The third term is 



T 4 = ,g 3 / dgS 4 (g) = 



7 



7 g 539 g 2 



2880 £ 3 768 S 2 1024 £ 



27101039 119119, f gE 

1 In — 

611440 2048 



The total result of the contour integral of the WKB expansion \sT = T + T 2 + T4, where 



T 



7 



1 



2880 E 3 24 E 15 g 



7 41 141 (7 15 2 
"768^ + 16 + ~ i? + Il6 + T i? 



539 16431 



-E- 



3043 



1024 E 128 
27101039 862501 



E 3 



1865 1155^ 



64 



61440 



128 



-£ 2 



5616029 . 



512 



-E 4 



16 



119119 285285 



2048 



256 



■£ 2 



255255 
128 



£ 4 In 



9E\ 
8 7 



The perturbative counterterms read 

P=|ln(27r)-ln [r (± - B 3 (E,g))] + B a (E,g) In 



7 



1 



79 




768 E 2 


f | + ( 


539 


-^+: 


1024 £ 4 


16 


"319081 


175325 


8192 


1024 



7_ 15 
16 + T 



225 
32~ 



£ 3 



128 



/1365 
\~64~ 

-£ 4 4 



1155 



2048 



£ d In — 



E z 



2880 £ 3 24 £ 



gE 



-E + E In 



256 



128 



E In 
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where the result of B 3 (E,g) from Eq. (|22|) has been used, as well as the asymptotic expansion 
In r(i + z) = z {ln(z) - 1} + A ln(27r) 



24z 2880 z 3 



Note that the logarithm \n(gE/8) is obtained as the combination of B 3 (E,g) ln(— g/8) ss E ln(— g/8) and from 
the asymptotic expansion of the expression — lnT(i — B 3 {E,g)) w E ln(— E). Finally, the instanton function 
A 3 (E,g) for the cubic potential is found to read 



A 3 (E,g)=T-P 



2 
155 



77 

32 



141 



9 



( 49415863 6724683 



V 122880 



1024 



E' 



E- 



9 



2 ^5911^ 



V 128 



11947 
32 



E 



5481929 
512 



0(g 4 ). 



(46) 



(47) 



This result for the instanton function is the basis for the calculation of higher-order corrections to the decay width 
of the resonance energies, as will be explained in the following. 



4 Modified quantization condition 

For g < 0, the interaction Hamiltonian ^/gq 3 = ±i\^/g\q 3 is 7>T-symmetric, and its spectrum is purely real. In 
view of Eqs. (|20|) and (f2Tj) . we then have the purely perturbative quantization condition 

B i {E,g)=n+\ & -jj 1 =0, g<0. (48) 

r{±- B 3 (E,g)) 

Let us now study how this condition needs to be modified in case instantons exist, i.e., for g > 0. We start from 
the unperturbed harmonic oscillator 

In this form, we naturally identify 

r(1 l p-, = det(« - E) (50) 
I (j - E) 

as the spectral determinant. The effect of the instanton is to add an infinitesimal imaginary part to the energy of 
the bound state, 

E ^ £ = n+\+ilm£ . (51) 

Indeed, the definition of resonance energies via complex scaling and the proof of their existence can be found in 
Ref. [7]. Expanding in the nonperturbatively small (in g) imaginary part of the energy, we obtain 

_(_!)n+l/2 n , Img (52) 



r{\-8) 

For the cubic oscillator, we approximate the right-hand side of Eq. ([52)) as 

lva£ « Ime^(g) ps -^—g- n - 1 ' 2 exp ( -t^- ) « A= (-) exp (-A 3 (i?, 5 )) . (53) 

wiV 7 !" V 15 3/ n!V87r \5/ 

We now generalize the left-hand side of Eq. (|52|) as 1/^(5 — B 3 (E,g)) and approximate the right-hand side of 
Eq. (|52|) by the right-hand side of Eq. (f53|) . and obtain 

B 3 (S,9) 

exp[-^a(E,a)], g>0, (54) 



r(|-B 3 (S, 5 )) v 7 ^ V 9 

which is our conjecture for the resonance energies of the cubic anharmonic oscillator in the "unstable" regime of 
positive coupling parameter g > 0. When the logarithm of both sides of Eq. (|54]) is taken, we recover the structure 
of the terms on the right-hand side of Eq. (f30|) . 

In Sec. [3l we have calculated the function A 3 {E,g) in the WKB limit g — > 0, with gE fixed. In order to derive 
the result (|46|) for the function A 3 (E,g), we had to expand the function B 3 (E,g) in the limit of large E, which 
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implies [according to Eq. (05])] the approximation — ln.T(| — B 3 (E,g)) m Eha.r(—E). A priori, this expansion is 
not applicable when E w n + A + 5, where 8 summarizes the perturbative and nonperturbative corrections (in g). 
We have calculated the instanton function A 3 (E,g) in an "unphysical" domain of large quantum numbers n (large 
E), "away" from the first few resonance energies. However, the result ([46]) is written in terms of an expansion 
whose terms actually decrease for typical values of E for the first few resonance energies, where E is close to a 
small half-integer. We can thus make an — in some sense fortunate — observation: While the expansion (|45|) for the 
logarithm of the r function is not applicable to the physical domain of E = n + | + S, the obtained expansion for 
the A 3 (E,g) function is applicable in both the "unphysical" WKB domain of gE fixed, and g — » 0, thus E — >• oo, 
as well as the physical domain of E « n+ 1, and g —> 0. We can thus use the modified quantization condition (|54|) 
in the regime of low principal quantum numbers, which is of prime interest. 

This unified result, which gives immediate and systematic access to the higher-order corrections, should be compared 
to other approaches (Refs. [13] [191 [22j [23] , which also strive to go beyond the simple leading-order results for the 
decay rates of the resonances. One approach [13] [T9] uses variational^ improved perturbation theory which leads 
to a numerical improvement over the simple summation of the first few perturbative terms in the strong-coupling 
domain, but does not lead to a systematic prescription for calculation of higher-order coefficients in the "resurgent" 
expansion. The other approach [22] [23] is based on a set of quantization conditions derived as variations of the 
Bohr-Sommerfeld condition. 



5 Results 



Based on the results presented in Eqs. ([22]) . ([46]) and ([54]) . it is easy to derive higher-order corrections to the 
imaginary part of the first few resonance energies of the cubic potential, by entering into the condition ([54]) with 
an ansatz 



o3n / o \ 

lme®(g) = —= g""" 1 / 2 cxp -— (l + a n g + b n g 2 + ...) , for g > , (55) 

and equating coefficients. The higher-order coefficients, for the lowest four states, read 

(3)/ s cx p( _ tI^) f 169 44507 , 86071851 , 189244716209 . m , 

Imei ; (o) = ^ \ \ g g 2 g 3 g 4 + O g 5 ) } , (56a) 

yy ' y/Wg \ 16 y 512 y 40960 y 2621440 y va ') ' v ; 

T (3) . , 8e- 2 /( 159 ) f 1 853 33349 2 395368511 , 1788829864593 4 ml 

Im£ ^ {9) = ~ \ 1 --l6-fl + -5iT^- ^0960- 5 2621440— 9 + ° {y } | ' (56b) 

Im4 3) ( g ) = - 32 % 2/ !; 5 2 9) (l - —9 + H^l, 2 - 1085785671 g 3 - 4272925639361 j 
V^ff 5/ I 16 512 y 40960 y 2621440 y yy ') ' 

(56c) 

(3)/ , 256e- 2 /( 159 ) f 3913 8807869 2 1571666861 „ 3214761534593 , m . 

Imei (g) = == < 1 g H g g g + O(q ) 

3 30Fg 7 /2 \ 16 y 512 y 40960 y 2621440 y yy ' 



(56d) 

For the coefficient of relative order g, a general result (for any n) has been indicated in Eq. (22) of Ref. [I]. 
For n = 2, Eq. (22) of Ref. [1] gives a correction term of —2041(7/16 instead of —2101(7/16 for the state with 
n = 2, and a correction term of —10543(7/48 instead of —3913(7/16 for the state with n = 3. With the help of the 
dispersion relation ([4]), these coefficients can be related to the corrections to the leading factorial asymptotics ([8]) 
of the perturbative coefficients of relative order l/K. We have checked our calculation against numerical values 
of the perturbative coefficients of up to 60th order, which can easily be determined on the basis of the relation 
B 3 (E,g) =n+l/2. 
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6 Conclusions 



The calculation of higher-order corrections to the decay widths of quantum states in unstable potentials is a 
challenging problem. Analytic approximations are useful in that regard, if they contain a sufficient number of 
correction terms which also permit an estimate of the truncation error. As shown in Sees. [4] and [5] the calculation 
of higher-order corrections to the decay widths of resonances of the cubic anharmonic oscillator, is based on two 
characteristic functions, namely, the 

perturbative function B 3 (E, g) = E + 0(g) defined in Eq. ([21]) and evaluated in Eq. ([22]) (57) 

and on the 

2 

instanton function A 3 (E,g) = — — \-Oig) defined in Eq. ((30)) and evaluated in Eq. ([46]). (58) 

15 g 

The first of these functions is given by an evaluation of a contour integral about the poles of the perturbative 
approximation to the logarithmic derivative of the wave function (see Sec. [2]). The second function, A${E,g), is 
determined by a contour integral about the cut of the WKB approximation of the wave function (see Sec. [3D. 

The modified Bohr-Sommerfeld quantization condition ((4]) then gives us access to the higher-order corrections to 
the imaginary part of the resonance energies. These are given, for the first few resonances of the cubic potential, 
in Sec. [5] The full solution of the modified quantization condition has the structure of a generalized nonanalytic 
expansion, or "resurgent" expansion, 



oo 



4 3) (<?) = £ 

J=0 



1 

"2 ( 2 
exp 



!V8^ \gj V 15 <? 



j 

£^H)E S Ki^> 9>0, (59) 



q , 

L=0 x a ' K=Q 



where 



£max = max(0, J — 1) . (60) 

The perturbative coefficients e^ K are equal to the Sq 3 ^ coefficients. The coefficients entering Eq. ([56]) are the 
3jq| coefficients. 

The coefficients 3p q i multiply the correction of relative order g to the one-instanton (J = 1) effect and are given 
in Eq. ([56]) for n = 0,1,2,3. They grow quite drastically in magnitude with the principal quantum number. In 
the notation of Eq. (|55p , we have for the state with n = 3 a coefficient of a 3 = —244.563. This means that 
the correction term of order g halves the total result for the decay width of the third excited state of the cubic 
potential already at a minuscule coupling parameter of g « 0.002. It is quite interesting to see that numerically 
approximations to such a fundamental physical quantity as the resonance energies of the cubic anharmonic oscillator 
require a considerable effort in their evaluation. Here, our aim was to carefully explain the calculational approach 
necessary in order to gain access to the higher-order corrections, which are numerically significant. 
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